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The Letter analyzes stability of spin precession currents in superfluid '^He-B when the precession 
angle is very close to 104°. In this limit a spin-precession vortex has a very large core, and a barrier 
blocking motion of these large-core vortices across the current streamlines (phase slip) disappears at 
precession-phase gradients much smaller than critical gradients estimated from the Landau criterion. 
Nevertheless, spin precession currents remain stable up to the Landau critical gradients, since in 
this case there is a barrier, which blocks the phase slip at very early stage of vortex-core nucleation. 
The Letter also predicts a second-order phase transition between a parity-symmetric and parity- 
nonsymmetric spin-precession vortex cores at the precession angle 126.5°. 
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The phenomenon of spin superfluidity was intensively 
studied experimentally and theoretically in 70s and 80s 
of the last century (see the review articlesi^i^ and ref- 
erences therein). Nowadays there is a revival of interest 
to the phenomenon of spin superfluidity in ^He-B^i^i^i^i^. 
Efforts to observe a similar phenomenon in '^He-A were 
also reported^. Meanwhile, there still remains some un- 
resolved controversy, which was inherited from the previ- 
ous rounds of studying this problem: what is the struc- 
ture of spin- precession vortex^*', which participates in the 
phase-slip process, and what determines stability of the 
spin-precession supercurrent in "^Hc-B. 

The concept of spin vortex as a phase-slip tool de- 
termining stability of the superfluid spin current was 
introduced in 1978^^. For the superfluid "^He-B the 
spin-precession vortex was considered in 1987— (see also 
Ref.^^ ). It was obtained that the core radius Tc of the vor- 
tex is on the order of the dipole length: rc S,d = c±/fl, 
where is the longitudinal-NMR frequency and c± is 
the velocity of transversal spin waves. Using the Landau 
criterion it was shown that the critical phase gradient is 
also determined by the inverse dipole length. The bar- 
rier for vortex growth in the phase-slip process vanishes 
at phase gradients of the order of the inverse core radius. 
So the threshold for vortex instability agrees with the 
critical gradient in the Landau criterion. This is usual in 
the superfluidity theory^. 

One year later Fomini^ suggested that the vor- 
tex core must be determined by another scale i^f — 
c±/ y/{u!p ~ uil)ul1l, where ujp and lol are the preces- 
sion and the Larmor frequencies. This was supported 
by Misirpashaev and VolovikiS' on the basis of the topo- 
logical analysis. Since l/^f. plays a role of the chemical 
potential for the precession moment and is directly con- 
nected with S^d, the question whether the core radius is 
determined by S^d or S^p, is similar to the question whether 
the core radius in the bose- liquid is determined by the liq- 
uid density, or by the chemical potential: Both answers 
are correct since the quantities are connected by thermo- 
dynamic relations. A real important difference was that 
according to Fomin if the precession angle (3 approaches 



to the critical value j3c = 1.82 rad (or 104°) the core ra- 
dius becomes Tc ~ ~ S,d/{P — Pc), i-e., by the large 
factor l/{(3 — f3c) differs from the estimation Tc ~ £,d done 
in Ref.J"^. So the latter is valid only far from the critical 
angle, where P — (3c ^ 1- 

Since no barrier impedes vortex expansion across a 
channel if the gradient is on the order of l/rc, the large 
core Tc ^ S,d/{P — Pc) at /3 — *■ /3c leads to the strange 
(from the point of view of the conventional superfluid- 
ity theory) conclusion: The instability with respect to 
vortex expansion occurs at the phase gradients ~ l/rc 
essentially less than the Landau critical gradient ~ 
obtained in Ref,-- for any j3 > (3c- The present Letter 
suggests resolution of this paradox. It demonstrates that 
at precession angles close to 104° at phase gradients less 
than the Landau critical gradient but larger than the in- 
verse core radius no barrier impedes phase slips at the 
stage of vortex motion across streamlines, but there is a 
barrier, which blocks phase slips on the very early stage 
of nucleation of the vortex core. So for these gradients 
stability of current states is determined not by vortices 
but by vortex-core nuclei. 

The analysis addresses also possible symmetries of the 
vortex core. It was expected that parity symmetry (its 
deflnition is given below) is always brokeiJ^. The present 
Letter presents numerical calculation demonstrating the 
second-order transition between a parity-symmetric and 
a parity-nonsymmetric vortex at the precession angle 
126.5°. In the past the first-order transition in cores of 
■^He-B mass vortices was detected in NMR experiments 
on rotating ^He-B^^. It was theoretically explained in 
RefsJ^ii^ in terms of the transition between the axisym- 
metric and non-axisymmetric cores. Later this theory 
was confirmed experimentally by direct observation of 
the non-axisymmetric core in one of the two vortices^^. 

The spin dynamics of superfiuid phases of "^He is de- 
scribed by the theory of Leggett and Takagi^S. Following 
Fomin^ we introduce the Euler angles a, (3, and 7 in the 
spin space of the '^He-B order parameter. The angle (3 is 
the precession angle, and a is the precession phase. The 
angle $ = a -f 7 characterizes the resultant rotation of 
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the order parameter in the laboratory frame, and in the 
hmit f3 ^ (no precession) becomes the angle of rotation 
about the z axis. The moments canonically conjugate to 
the angles a, P, and <I> are respectively: P — Mz — M^, 
Mg, and M^, where is the z component of the mag- 
netization M in the laboratory coordinate frame, is 
the projection of M on the ^ axis of the rotating coordi- 
nate frame, and Mp is the projection of M on the axis 
perpendicular to the axes z and ^. 

For phenomena observed experimentally only one de- 
gree of freedom is essential, which is connected with the 
conjugate pair "precession phase a-precession moment 
P" . The Hamilton equations for the precession mode 
are: 



da 



SF 



dt ^SP' 



dP 
'dt 



Sa 



(1) 



Since the degree of freedom connected with the conjugate 
pair is not active, the angle $ is determined from 

minimization of the energy: = 0. The free energy 

F = Fz + -Fv + V includes the Zeeman energy Fz — 
—M ■ H = —MHu = — xw£/7^, where H = Hz is an 
external constant magnetic field, the gradient energy (we 
assume that the spin current is normal to the magnetic 
field Hz), 



F^ 



where 



A{u\ 



A{u) = -^{\~uf + 1-V? 



, (2) 



(3) 



and the dipole energy V = xc±v{u, ^)/7 Cd' where 
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w(u, $) 
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(1 -I- cos$)u -I- cos$ 



1 



(4) 



Here x is the magnetic susceptibility, 7 is the gyromag- 
netic ratio, ljl = iH is the Larmor frequency, u = cos /3, 
and the ratio c\\/c±_ of velocities of longitudinal and 

transversal spin waves will be chosen to be ^/aJ?}^. In 
the state of stationary precession the precession angular 
velocity is constant: da/dt = —ujp. This state corre- 
sponds to the extremum of the Gibbs thermodynamic 
potential, which is obtained from the free energy with 
the Legendre transformation G — F + ojpP/^. Thus the 
precession frequency ojp plays the role of the "chemical 
potential" conjugate to the precession moment density 
P. The distribution of the parameters u — cos (3 and 
$ is determined from the two Euler-Lagrange equations 
5G/6u ^ and 5G/6^ = 0. 

For uniform precession minimization with respect to u 
yields the relation 
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1 dv 



r 



(5) 



and minimization with respect to $ (only the dipole en- 
ergy depends on $) gives the equation 



(1 -f cos $)u -f cos * — ^ 



(l + u)sin$ = 0. (6) 



Solution of this equation yields 
1/2 -u 



cos$ 



, v{u, $) = 



l + u 

for /3 < 104° [u > -1/4) and 



cos $ = 1, v{u,<^) = vq{u) ^ — [^ + u 



(7) 



(8) 



for /3 > 104° {u < -1/4). 

The spin-precession vortex state is non-uniform, and 
the gradient energy becomes essential. For an axially 
symmetric vortex with 27r circulation of the precession 
phase a (Va = 1/r) the Euler-Lagrange equations are 
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A-u 
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u / du 

(1-^2)2 
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(1 + cos <i>)u + cos $ - i 



(l + cos$) =0, (9) 



1 d$ 

r dr 



d^ 

dr^ 



1 



(1 -I- cos $)m + cos $ - - 



{l + u) sin$ = 0.(10) 



There are two types of vortices corresponding to two 
types of symmetry. The solution with $ = is parity- 
symmetric, while in the structure with $ 7^ symmetry 
with respect to parity transformation (f> — s- — $ is broken. 
At the periphery of the vortex core, where $ ^ 1 and 
u -I- 1/4 < 1, the solution of Eqs. ^ and is 
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$ = C 



64r2 



3$2 



r 



(11) 



where the constant C is zero for the symmetric vortex 
and is of the order of unity for the non-symmetric vor- 
tex. The ratio between Fomin's and the dipole length 
is determined by the value of u at infinity: = 
4Cf-\/|uoo + 1/4|/15. Equation (fTT|) demonstrates that 
in the hmit S^d £,f (I" + 1/4| <C 1) the first two terms 
in the expression for u can be neglected everywhere ex- 
cept for very large distances r of the order or larger than 
ln(^F/Cci)- It is the approximation of Fomin, who used 
the relation Eq. ([7]) between $ and u = cos (3 obtained 
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for the uniform state for u > —1/4. This allows to re- 
duce two coupled equations (O and pU)) to a single one 
after exclusion from them the terms oc l/Q^- The result- 
ing equation does not contain the dipolc length explicitly, 
but the length ^i?, which determines the core size, cer- 
tainly depends on it. 

In order to find whether the symmetric $ = vortex 
can be realized we studied stability of the solution of 
Eq. ([9]) for <i> = 0. This requires an analysis of a linear 
equation for <f>, 



1 d$ 

r dr 



d^ 

dr"^ 



2(m + 1/4)(1 + m) 
5d 



$ = e$, (12) 



which is analogous to the 2D Schrodinger equation for 
a particle in a potential well. Here u is determined 
from Eq. ([HI) at $ = 0. The <I> = vortex becomes 
unstable if Eq. has a solution with negative en- 

ergy e < 0. It is known from quantum mechanics that 
the 2D potential well has always a bound (localized) 
state^?^ But since for u < —1/4 (the necessary condition 
for stable spin-precession current) the "potential energy" 
2(u -I- 1/4) (1 -I- u)/5Q in Eq. is positive at infinity, 
the existence of a localized state leads to the condition 
e+2(uoo + l/4)(l+Uoo)/5^^ < 0, which does not rules out 
that e > 0, i.e., the solution $ = can be stable (except 
for the limits Uoo ^ — 1/4 and Uoo — > —1). The phase 
transition between two types of vortices is determined by 
the condition e = —2{uoo + 1/4)(1 -I- Uoo)/5Cd- Studying 
Eq. p2)) numerically we have found that e reaches its 
critical value at Uqo = —0.594, which corresponds to the 
precession angle f3 — arccos(uoo) — 126.5°. Using the so- 
lution for $ in the critical point for estimation of the en- 
ergy contribution of the terms of order we have found 
that this contribution is positiv622i. Therefore the phase 
transition between the symmetric and the nonsymmetric 
vortices is of the second order. The numerical solutions 
of Eq. (IH) (with $ = 0) and Eq. are shown in Fig. 1. 
We analyzed the spin-precession vortex in the absence 




FIG. 1: FIG. 1: Vortex structure at the transition from 
parity-symmetric to parity-nonsymmetric vortices. The plot 
of <& (dashed line) is normalized to the value "I'(O) = 1. 



of spin current at infinity. If an uniform spin current 
is present, motion of vortices across current streamlines 
leads to precession-phase slips suppressing the current. 
If the precession-phase gradient connected with the cur- 
rent exceeds the inverse core radius, there is no barrier 
impeding this process. The analysis of the Landau crite- 
rion showedi^ that the spin-precession current is stable 
up to the critical gradient on the order of the inverse 
dipole length l/^d- In the limit u — > —1/4 the latter can 
essentially exceed the inverse core radius y^\u + l/4|/^c;. 
But how can currents with gradients between and 
\/\u + ^T^\/£.d be stabihzed? The answer to this ques- 
tion is that the barrier blocking phase slips is present 
at the very early stage of nucleation of vortex rings in 
3D systems (or vortex pairs in 2D systems). Vortex nu- 
cleation starts from a slight localized depression of the 
superfluid density (determined by A{u) in our case). In 
order to analyze this " protonucleus" of the vortex core in 
the spin-current state one should consider a new Gibbs 
thermodynamic potential G — G — xc\j'^a/j^: 



G 



2{l-u^) 



v{u) 

^d 



(13) 



where j is a Lagrange multiplier. The nucleus, which is 
related with a peak of a barrier, corresponds to a saddle 
point in the functional space. So its structure should be 
found from solution of the Euler-Lagrange equations for 
the introduced Gibbs potential. The first step is to vary 
the Gibbs potential with respect to a. Let us restrict 
ourselves with a ID problem, when the distribution in 
the nucleus depends only on one coordinate x. Then 
the distribution of Va is given by Va = j/A{u), where 
the reduced spin-precession current j = A(Moo)Vao is 
determined by the gradient Vag far from the nucleus 
center. Expanding with respect to small deviation g — 
u — Uoo from the equilibrium value of u at infinity one 
obtains 
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G 



] 



2Aiu) 2(1 -«2) 



(V.9)2 



2(1 -u 
" 2 du^ 



v{Uc 



2A{Uao) 
w(Uoc 



v{u) 

1 

w 



2A(woo) Q 



6 2 



du^ 



-(14) 



where we took into account that d^v{u)/du^ = 0. The 
linear in g terms must vanish at the stationary current 
state. The term quadratic in g determines the stability 
of the current state: it vanishes at the Landau critical 
current 



Jc 



1 d^viuoo) ( d'^[A{uoo) 



du^ 



dv? 



(15) 



which was derived in Ref.'^^. Considering the case of the 
current close to the critical value and using the Taylor 
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expansion of A ^ (u) around u 



-1/4 one obtains 



G 



1572 



where 



a = 0.239(j2-j2), 6 = 0.577j2. 



(16) 



(17) 



The Euler-Lagrange equation for this Gibbs potential, 
— Ag = ag — bg^ /2 = 0, determines the distribution of g: 



5 = 30 1 - tanh 



(18) 



where g^ = "ia/b = 1.24(j^ — j"^)/ i1 is the value of g in 
the nucleus center and rp = = 4.1^d/-\/jc ~ P is 

the nucleus size. The energy of the nucleus, 
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3a/b 



6g3 64xci a5/2 



0.2US 



Xcl [Jl-Pf^ 



J2 



,(19) 



determines the barrier for the process of the vortex core 
nucleation. Here 5* is the cross-section area of the chan- 
nel. Since in the limit j — *■ jc the nucleus size Tp is 
divergent our ID description is always valid close enough 
to the critical point, where rp 3> ^/S. When Vp be- 
comes smaller than the transverse size of the channel, 
one should consider the 3D or 2D (in the case of a thin 
layer) nucleus. The first stage of this problem is to 
find the distribution of Va from the continuity equation 
V[(7Va] — 0. Its solution demonstrates that outside the 
nucleus the distribution of Va is the same as around the 
vortex ring (3D case) or the vortex dipole (2D case). In 
particular, in the 2D case 



Va — Vag — 



g{ri)rl dri 



Vao 2f(f-Vao) 



(20) 



In contrast to the ID case, the relation between Va and 
g is not local, so the following variation of the Gibbs po- 
tential with respect to g leads to an integro-differential 
equation. However on the basis of our solution of the ID 
problem and using the scaling arguments one may con- 
clude that the nucleus size can be roughly estimated from 
the expression Eq. (fT9|) for the ID case with replacing 5* 
by Tp or by rpd for the 3D and the 2D case respectively 
(d is the thickness of the 2D layer). 

In our analysis we assumed that the value of the pre- 
cession angle f3oo = arccos Mqo far from centers of a vortex 
or a nucleus was fixed and exceeded 104° (uoo < ~l/4). 
Otherwise (uao > —1/4) the dipole energy vanishes, 
and without dipole energy no stable current is possible. 
Meanwhile, Fomin^^ suggested that the spin-current can 



be stable even if m > —1/4 and the Landau criterion is 
violated. He argued that emission of spin waves, which 
comes into play after exceeding the Landau critical gra- 
dient, is not essential in the experimental conditions (see 
also a similar conclusion after Eq. (2.39) in the review 
by BunkoA*^). This argument is conceptually inconsis- 
tent. If the experimentalists observed "dissipationless" 
spin transport simply because dissipation was weak, it 
would be ballistic rather than superfluid transport. The 
essence of the phenomenon of superfluidity is not the ab- 
sence of sources of dissipation, but ineffectiveness of this 
sources following from energetic and topological consid- 
erations. The Landau criterion is an absolutely necessary 
condition for superfluidity. Fortunately for the superflu- 
idity scenario in '^He-B, Fomin's estimation of the role 
of dissipation by spin-wave emission triggered by viola- 
tion of the Landau criterion is not conclusive. He found 
that this dissipation is weak compared to dissipation by 
spin diffusion. But this is an argument in favor of impor- 
tance rather than unimportance of the Landau criterion. 
Indeed, spin-diffusion, whatever high the diffusion coeffl- 
cient could be, is ineffective in the subcritical regime, in 
which the gradient of the "chemical potential" wp is ab- 
sent. On the other hand, in the supercritical regime the 
"chemical potential" is not constant anymore and this 
triggers the strong spin-diffusion mechanism of dissipa- 
tion. 

It is worthwhile to remind that the Landau criterion 
is a necessary but not sufficient condition for transport 
without dissipation. At the Landau critical gradient (cur- 
rent) the current state ceases to be metastable because 
the barrier leading to metastability vanishes. But it is 
well known that in reality current dissipation via phase 
slips is possible even in the presence of barriers (due to 
thermal fluctuations or quantum tunneling). Therefore, 
the present work addresses only "ideal" critical currents 
(the upper bound for them) leaving "practical" critical 
currents beyond the scope of the analysis. 

In summary, the Letter analyzed stability of spin pre- 
cession currents in superfluid ^He-B when the preces- 
sion angle /3 is very close to /3c = 104° and the spin- 
precession vortex has a core of the radius ^ £,d/ V ~ Pc 
much larger than the dipole length S^d- Though a barrier 
for motion of these large-core vortices across the current 
streamlines disappears at rather small precession-phase 
gradients ^ \/ (i — l^c/id, spin- precession currents remain 
stable up to much large gradients ^ 1/^d, which were 
estimated from the Landau criterior^i^. Stability of cur- 
rents in this case is provided by barriers at very early 
stage of vortex-core nucleation. It was also demonstrated 
that at the precession angle 126.5° there is a second-order 
phase transition between a parity-symmetric and parity- 
nonsymmetric spin-precession vortex cores. 

The author thanks I. A. Fomin and G. E. Volovik for 
critical comments, which influenced conclusions of the 
present Letter. 



5 



^ E. B. Sonin, Superflows and superfluidity, Usp. Fiz. Nauk 
137, 267 (1982) [Sov. Phys.-Usp. 25, 409 (1982)]. 

^ I. A. Fomin, Spin currents in superfluid '^He, Physica B 
169, 153 (1991). 

^ Yu. M. Bunkov, Spin supercurrent and novel properties of 
NMR in ^He, in Progress of Low Temperature Physics, 
vol. 14, edited by W. P. Halperin, Elsevier, Amsterdam, 
(1995), p. 68. 

* G. E. Volovik, Twenty years of magnon Bose con- 
densation and spin current superfluidity in '^He-B, 
,cond- mat/0701180 

^ Yu. M. Bunkov and G. E. Volovik, Magnon condensation 
into a Q baU in ^He-B, Phys. Rev. Lett. 98, 265302 (2007). 

® Yu. M. Bunkov, Spin supercurrent, J. Magn. Magn. Mater. 
310, 1476 (2007). 

Yu. M. Bunkov and G. E. Volovik, Bose-Einstein conden- 
sation of magnons in superfluid '^He, J. Low Temp. Phys. 
150, 609 (2008). 

* Yu. M. Bunkov and G. E. Volovik, Spin vortex in magnon 
BEG of superfluid ^He-B, Physica C 468, 135 (2008). 

^ T. Sato, T. Kunimatsu, K. Izumina, A. Matsubara, 
M. Kubota, T. Mizusaki, and Yu. M. Bunkov, Observa- 
tion of coherent precession of magnetization in superfluid 
^He A-phase, arXiv: 0804.2994. 
^° As well as in Refs. IT^Il^ we prefer to tell about spin- 
precession rather than spin vortex since strictly speaking 
in ^He-B one deals not with currents of spin but with cur- 
rents of precession moment (see further in the Letter). 
E. B. Sonin, Analogs of superfluid currents for spins and 
electron-hole pairs, Zh. Eksp. Teor. Fiz. 74, 2097 (1978) 
[Sov. Phys.-JETP 47, 1091 (1978)]. 

E. B. Sonin, Superfluid transport of precession in '^He-B, 
Pis'ma Zh. Eksp. Teor. Fiz. 45, 586 (1987) [JETP Lett. 



45, 747 (1987)]. 

E. B. Sonin, Magnetic superfluidity in ^He, Zh. Eksp. Teor. 
Fiz. 94, 100 (1988) [Sov. Phys.-JETP 67, 1791 (1988)]. 
I. A. Fomin, Steady-state spin current in '^He-B, Zh. Eksp. 
Teor. Fiz. 94, 112 (1988) [Sov. Phys.-JETP 67, 1148 
(1988)]. 

T. Sh. Misirpashaev and G. E. Volovik, Topology of coher- 
ent precession in superfluid ^He-B, Zh. Eksp. Teor. Fiz. 
102, 1197 (1992) [Sov. Phys.-JETP 75, 650 (1992)]. 
^« O. T. Ikkala, G. E. Volovik, P. J. Hakonen, Yu. M. Bunkov, 
S.T. Islander, and G. A. Kharadze, NMR on rotating su- 
perfluid ^He-B, Pis'ma Zh. Eksp. Teor. Fiz. 35, 338 (1982) 
[JETP Lett. 35, 416 (1982)]. 

E. V. Thuneberg, Identification of vortices in superfluid 
^He-B, Phys. Rev. Lett. 56, 359 (1986). 
^® M. M. Salomaa and G. E. Volovik, Vortices with sponta- 
neously broken axisymmetry in '^He-B, Phys. Rev. Lett. 
56, 363 (1986). 

Y. Kondo, J.S. Korhonen, M. Krusius, V.V. Dmitriev, 
Yu. M. Mukharsky, E. B. Sonin, and G. E. Volovik, Direct 
observation of the nonaxisymmetric vortex in superfluid 
^He-B, Phys. Rev. Lett. 67, 81 (1991). 
A. J. Leggett and S. Takagi, Orientational dynamics of su- 
perfluid ■^He: A "two-fluid" model. I. Spin dynamics with 
relaxation, Ann. Phys. (N. Y.) 106, 79 (1977). 
LA. Fomin, Critical superfluid spin current in ^He-B, 
Pis'ma Zh. Eksp. Teor. Fiz. 45, 106 (1987) [JETP Lett. 
45, 135 (1987)]. 

L.D. Landau and E. M. Lifshitz, Quantum mechanics 
(Pergamon Press, Oxford, 1977), Sec. 45. 
Estimating the terms oc one must take into account 
corrections of the order to u. 



